Appli'-d  hese'.r'V. 


■ 


DETERMINATION  OF  TEMPERATURE,  VELOCITY,  AND  CONCENFRATI  ?N  PROFJ 


IN  THE  MI*  INC  I A  i  £R  BETWEEN  A  ROCKET  EXikt  F  1  I  A’.j 


"HE  SURj  JliNLlNG  SUPH  NIC  AIR  STREAK 


III.  J  PROCEDURE  FOR  COMPUTATIONS  AND 


PRELIMINARY  RESULTS  FOR  A  TYPICAL  CASE 


J.  Vasili u 
23  March  I960 


i 


Distribution  of  tni> 

U  villalUd. 


u.t 


1  cocumcat  .  uc  b<(  n.  opprovod 

for  public  iclooso  and  rale;  Its 

d’Tf  .button  is  unlimited _ ENGINEERING  DEPARTMEN1 


GO 


Is 


l  D*** 


CONVAIR  (ASTRONAUTICS)  DIVISION 
GENERAL  DYNAMICS  CORPORATION 


Reproduced  by  the 

CLEARINGHOUSE 
for  Foaoraf  Scientific  &  Technical 
Information  Springfield  Va.  22151 


C0NVMO. 

ASTpn* 


APR  19  1960 


BEST 

AVAILABLE  COPY 


EnR-AJM &jf  l a. 

Apjij  1  •..••I  He.-  ea.'  C  i 


DP,'j  SEMINATION  OF  TEMPERATURE,  VELOCITY,  AND  CONCENTRATION 
PROFILES  IN  THE  MIXING  LAYER  BETWEEN  A  ROCKET  EXHAUST 
JET  AND  TEE  SITRROUNDING  SUPERSONIC  AIR  STREAM 

PART  III.  PROCEDURE  FOR  COMPUTATIONS  AND 
PRELIMINARY  RESULTS  FOR  A  TYPICAL  CASE 


J.  Vasiliu 

23  i’.’oO 


iiUirlbution  of  thll  d,oouffllQi 
K  unlimited. 


I 


K.» 'INFER]'!'.'  i  JVPAi'r'MvJt 


in  iiar*.  untie).  Oor.v»fi  y  a 
b n  i  ili  pnvt  mtirf  ConUT 
.irirA  Orlny  ’.2~  '? 


:-3 


"i 


1R  (ASThONAR 
'Ll  !>)N/\K1CI 


■  )nl  Ob 


CONTENTS 


Pago 

ABSTRACT  . 

I.  DETERMINATION  OF  THE  INITIAL  AND  BOUNDARY  CONDITIONS . i 


A.  Initial  Conditions  . 

B.  Boundary  Conditions  . 

1.  Upper  Boundary  .  5 

2.  Lower  Boundary  . 

C.  Relations  for  the  Pressure  and  Pressure  Gradient  ....  7 

II.  DATA  REQUIRED  IN  THE  CALCULATIONS . 8 

A.  Thermodynamic  Functions  ...  .  8 


B.  Reference  Values  and  Dimensionless  Parameters  .  11 


C.  Input  Data . 1 4 

III.  COMPUTER  PROGRAM . lc 

IV.  TYPICAL  CASE  AND  PRELIMINARY  RESULTS . io 

REFERENCES . 26 

NOMENCLATURE . £8 

APPENDIX . 30 


ILLUSTRATIONS 


Figure  Title  Page 

1.  Schematic  Diagram  of  the  Mixing  Layer  . 


2.  Velocity  Profiles  During  the  Initial  Stages  of 

Development  of  the  Mixing  Layer . 21 

3.  Dimensionless  Velocity  versus  the  Dimensionless 

Coordinate  y* . 21 

4.  Dimensionless  Temperature  Frofiles . '22 

5.  Dimensionless  Temperature  versus  the  Dimensionless 

Coordinate  y* . 23 

6.  Concentration  Profiles  During  the  Initial  Stages  of 

Development  of  the  Mixing  Layer . 

7.  Computed  Thickness  of  the  Mixing  Layer  . 

8.  Experimental  Data  on  the  Value  of  0  . 


i 


v r  DETERMINATION  OF  TEMPERATURE,  VELOCITY,  AND  CONCENTRATION 

PROFILES  IN  THE  MIXING  LAYER  BETWEEN  A  ROCKET  EXHAUST 
JET  AND  THE  SURROUNDING  SUPERSONIC  AIR  STREAM 

PART  III.  PROCEDURE  FOR  COMPUTATIONS  AND 
PRELIMINARY  RESULTS  FOR  A  TYPICAL  CASE 

J.  Vasiliu 

ABSTRACT 

In  Section  I  of  this  Part,  the  method  of  obtaining  the 
initial  ar d  boundary  conditions  for  available  inviscid  Jet 
flow  solutions  is  discussed.  Simple  relations  for  the  determi¬ 
nation  of  the  temperatures,  velocities,  Mach  numbers,  pres¬ 
sures,  and  pressure  gradients  are  derived.  In  Section  II  the 
thermodynamic  functions  required  in  the  calculations  are  dis¬ 
cussed  in  detail  and  listed  in  Table  I.  A  program  written  in 
Fortran  II  that  makes  it  possible  to  carry  out  the  solution 
***■  of  the  problem  on  an  IBM-704  computer  is  described  briefly. 

Finally  one  typical  case  is  considered  and  plots  of  velocity, 
temperature,  and  concentration  profiles  are  given  for  the 
initial  stages  of  development  of  the  mixing  layer. 

I.  DETERMINATION  OF  THE  BOUNDARY  CONDITIONS 

It  was  pointed  out  in  Section  IV  of  Part  II  that  values  of  the 
velocities  and  temperatures  at  the  upper  and  lower  boundaries  may  be 
obtained  from  the  inviscid  flow  solution  of  the  Jet.  Such  a  solution 
was  obtained1  by  the  method  of  characteristics  for  the  case  of  an 
axially  symmetric  Jet  of  gas,  with  constant  ratio  of  specific  heats, 
issuing  into  a  supersonic  air  stream.  For  lack  of  other  data,  the  results 
of  this  analysis  will  be  used  to  obtain  the  boundary  conditions  for  the 
two-dimensional  mixing  layer  considered  here. 

A.  Initial  Conditions 

The  conditions  at  x  =  0  require  a  knowledge  of  the  velocity  and 
temperature  of  the  Jet  and  the  air  stream  at  a  section  normal  to  the 


1 


boundary.  The  flow  parameters  available  from  the  inviscid  flow  solution1 
are  indicated  in  Figure  1.  It  is  noted  that  the  Mach  number  of  the  jet  at 
a  section  normal  to  the  boundary,  denoted  by  M  ,  is  different  from  the 
Mach  number  at  a  section  normal  to  the  exit,  denoted  by  M.  .  ,  the  former 
being  larger  in  view  of  the  expansion  taking  place  at  the  periphery  of 
the  nozzle.  The  initial,  conditions  will  be  based  on  M  •  Using  the  total 
temperature  of  the  jet,  Tip  ,  and  the  constant  ratio  of  specific  heats, 

TT^,  one  has 

T 

T, 


r 1t-i 

1  +  M2 


(1) 


The  velocity  is  expressed  in  terms  of  the  Mach  number  and  temperature 
from  the  relation 


U 


M 


r..gR  t 

jtu  m  e 
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(2) 


Where  R  is  the  gas  constant  for  the  mixture  at  the  exit  conditions  given 
by  *e 


R  =  — 
m  W 
e  ra 


(3) 


with  W  ,  the  molecular  weight  of  the  mixture  at  the  exit,  given  by 


Wm  -  ZWA 

G  i 

Substituting  Eq.  (j)  into  (2),  and  using  consistent  units,  one  findr 


U  =  222.97  M 
e  e 


Jrj£i 

M  W 
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(ft/ sec) 


(4) 


For  the  air-stream  side  one  has,  from  the  energy  equation,  using  the 
total  temperature  of  the  stream  T^ 
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and. 
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The  velocity  con  be  obtained,  ac  before,  from  a  relation  analogous  to 
Eq.  (U) 


X  T 


Ua  =  222.97  M  h~  (ft/ r.ec) 
'  m 

n 


(6) 


v/here  U  ,  the  molecular  v/eight  of  the  air  stream,  ir.  Livon  by 
L  a 


U  = 
m 


‘  I 


,  (ci) 


(7) 
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The  concentration  of  species  at  the  nozzle  exit  can  be  obtained  from  an 
analysis  of  the  problem  of  chemical  reactions  during  the  expansion  of 
premixed  gases  through  the  nozzle.2  It  is  indicated  here  that  one  may 
consider  the  chemical  reactions  in  the  flov  to  be  in  local  equilibrium 
up  to  a  certain  point  of  the  expansion,  and  to  be  frozen  at  those  values 
thereafter,  up  to  the  nozzle  exit  area  and  further  downstream .  Such  a 
treatment  yields  values  of  the  specific  impulse  v/hich  are  in  good  agree¬ 
ment  v/ith  observations  on  actual  nozzles.2  The  three  additional  species, 
namely,  li,  0,  and  Oil,  listed  by  lieu,  and  Boynton* ,  are  neglected  in  this 
uork.  The  mole  fractions  of  the  constituents  used  in  the  calculations 
are  given  in  Section  III.  For  the  air  stream,  the  standard  argon-free 
air  composition  vill  be  used. 


B.  Boundary  Conditions 

The  method  of  characteristics1  gives  the  Mach  number  of  the  flov/  at 
either  side  of  the  Jet  boundary  in  terns  of  the  dimensionless  axial 


k 
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distance,  X/rc,  and  the  dimensionless  radial  distance,  r/rc.  It  is  easily 
Eeen  from  Figure  l.b  that  the  dimensionleGG  coordinate  x  used  in  this 
work  can  be  related  to  the  above  coordinates  (after  linearizing  each 
small  segment  on  the  boundary)  by  means  of  the  relation 


where  "i"  denotes  the  points  at  which  the  values. of  M.  and  M  are  known 

J  g 

(see  Table  II )  •  For  calculations  on  the  computer,  the  values  of  Mj  and 
M  may  be  listed  as  a  function  of  (x).,  or  thece  Hach  numbers  can  be 
expressed  as  polynomials  in  x.  It  has  been  found  that  the  maximum  error 
incurred  in  using  a  4th  deg  ?e  polynomial  approximation  does  not  exceed 
0.30  percent.  Values  of  the  coefficients  for  a  specific  case,  are  given 
in  the  section  on  input  data. 


1 .  Upper  Boundary 

In  view  of  the  above,  the  Mach  number  on  the  Jet  side  will  be 
given  by  the  relation 

Mj  -  aQ  +  OjX  +  a2x2  +  a3x*  +  a4x4  (9) 

The  temperature  at  this  boundary  can  then  be  expressed  in  terms  of  this 
Mach  number  as 


j.. - y  i 

or,  in  dimenBionless  form  and  using  the  notation  of  Part  II  for  the  grid 
point g  at  the  boundary, 

TT  /To 

*  A  A  A  A  u 

T  S  T  H  — «*•  SS  — — w .  ( 11 ) 
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The  velocity  is  obtained  by  a  relation  analogous  to  Eq.  (4),  i  .c . 


U.  =  222.97  H 
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J  >  V 
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which  can  be  rewritten  in  the  following  dimensionless  fonn 
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It  is  noted  that  in  Eq.  (l2)  and  (13)  the  molecular  weight  of  the  mixture 
at  the  upper  boundary  iB  the  Banc  as  that  at  the  nozzle  exit.  This  is 
consistent  with  the  remaininc  boundaiy  condition  i  .e . 


C 


i 

1J+1,  l 


(i4) 


which  iB  based  on  the  assumption  that  the  composition  of  the  Jet,  down¬ 
stream  from  the  nozzle,  remains  frozen  at  the  conditions  at  the  exit . 
Equation  (ll),  (13)  and  (l4)  constitute  the  bound-  ry  conditions  for  the 
uppe  •  boundary  of  the  raixirc  layer. 

2.  Lower  Boundary 

For  the  air- stream  side,  the  Mach  number  will  be  obtained 
from  the  relation 

M  =  P  +  P  x  +  P^x2  +  P  x3  +  P  x4  (ll;) 

B  CO  1  U  3  4 


Based  on  thin  Mach  number,  the  temperature  at  this  boundary  iB  Given  by 
a  relation  annloGOUs  to  Eq.  (5) 


1  + 


r 


(.16) 


or,  in  dimensionless  form,  and  in  terms  of  the  notation  for  a  Grid  point 
on  the  lower  boundary, 


r -1 


t  n  + 

t  L  r -1 

0  1  + 


(17) 


The  velocity  can  be  determined,  from  a  relation  analoGoue  to  Eq.  (6): 
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or. 


U  M  fr  T  A 

u_  =  —  =  PPP.97  \/—  — *—■ 

uO ,1  U  ^  ''  IJ  v  V, 

’  o  on 


a 

where  V?  is  Riven  by  I!q.  (7)*  The  additional  condition  on  the  concen¬ 
trations*1^ 


°0 ,1  Vt  a 


(PO) 


Equation  (17),  ( 19 )  and  (PO)  constitute  the  boundary  conditions  for  the 
lower  boundary  of  the  nixing  layer. 


C .  Relations  for  the  Pressure  and  Pressure  Gradient 

The  inviseid  flow  solution1  yields  the  stagnation  pressure  at  both 
sides  of  the  boundary.  Values  of  the  static  pressure,  required  in  this 
work,  can  be  obtained  fron  this  pressure.  Since  the  static  pressure 
Must  bo  balanced  at  the  jet  boundary,  one  r,ioy  use  the  flow  parameters 
at  either  side  in  carrying  out  the  calculations .  Selecting,  arbitrarily, 
the  stream  side  stagnation  pressure,  Pp  ,  and  Mach  number,  Mo,  the  static 
pressure,  p  ,  is  given  by 


P 


r  /(r.-i) 


(pi) 


or,  in  dimensionless  form, 


rjir  -i) 


r„-i  \ 


(pp) 


The  dimensionlecs  pressure  gradient  in  obtained  by  d i f f erent i at i ng  the 
above  e:cprcssion 


r  M  p 

yJ  U 


dx 


lc  ) 
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Fron  Eq.  (l^) 


+  2P2x 


+  5P 


3“ 


(?Jt) 


Sub  st  itut  ing  Eq . 


(2>v)  into  Eq.  (25), 


the  precDure  gradient  is  found  to  bo 


P 


x 


+  2P,  x  +  3P_xa  +  Mi . 

ti  J  »i 


(25) 


Tliis  equation  is  to  be  calculated  at  the  point  x  =  x  as  required  in 
Part  II .  The  pressure  at  the  nozzle  exit  can  be  obtained  from  the  :  tar¬ 
nation  pressure  of  the  jet .  Thus 


P 


c 


7r 


Jt 


-1 


(26) 


II  .  DATA  REQUIRED  III  THE  CAJCULATI01  IS 


A.  Thermodynamic  Functions 


In  Part  II  a  function  f^(T),  Eq.  (36),  has  been  UGed  in  the  equation 
for  the  dimensionless  specific  heat  of  each  species  to  denote  the  temper¬ 
ature  dependent  molar  specific  heat  of  the  can.  This  function  can  be 
expressed  by  means  of  approximating  polynomials  such  as  those  civ<  n  by 
Sueiccrt  and  Beardsley.  The  equations  arc  of  two  General  forms,  and 
they  are  valid  within  a  Given  rance 


f 


i 


c 


or 


ft(T)  =  a 


+ 


c 

T 


where  T  is  the  temperature  in  8R.  The  equations  eorrespondiuG  to  each 
species  are  listed  in  column  8  of  Table  I .  Their  accuracy  was  checked 
against  the  values  Given  in  the  Gas  Tables  for  the  ranee  of  500-3000°R. 
The  maximum  error  found  is  0  /)0  percent  for  the  CO,, -curve,  0 .80  percent 
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for  the  0  -  and.  II2- curve,  1.0  percent  for  the  CO-cu^ve,  1.5  percent  for 
the  K^,  and  2.0  percent  for  the  II20  curve.  It  in  noted  that,  in  Table  I, 

T  in  the  dimensionless  temperature . 

For  temperatures  from  500°R  <io\m  to  ubout  250°R,  the  specific  heat 

nay  be  assumed  constant.  Values  taken  from  the  Gas  Tables  arc  also  listed 

in  column  8.  For  0a,  !L,0,  CO,  IE,  this  assumption  leads  to  negligible 

error.  For  II  and  CO  the  constant  values  used  were  selected  arbitrarily 
2  2 

so  as  to  give  small  errors  (a  maximum  of  2.5  percent)  above  500°R.  Below 
this  temperature,  the  error  for  these  two  gases  is  appreciable,  and  it 
nay  become  necessary  to  use  approximating  equations. 

The  function  f*(T)  tabulated  in  column  9  is  the  derivative  of  f^(T) 
with  respect  to  the  dimensionless  temperature;  it  appeal's  in  some  terns 
of  the  energy  equation. 

A  method  of  calculating  the  dimensionless  enthalpy  of  each  species 
is  discucscd  Jr.  Port  II.  The  equation  obtained  there  [Kq.  Oil)],  is 


i  po  o  L  i  J 


Here  G^(t)  is  the  integral  of  f^(T),  and  is  given  in  colui.in  10.  Below 

500°R  this  term  in  set  equal  to  zero  as  noted  in  Part  II .  Hie  constant 

h  was  computed  from  the  relation  [Eq.  (4o)  of  Part  II ] 

‘  i 


b. 


000 

=  J  f  (T)dT  -  G 


(500) 


T  >  500° R 


For  temperatures  below  500°R  h^  is  given  by  Eq.  (42),  i  .c . 


h  =  f  (T)T 
i 


T  <  500°  R 


Tlic  above  values  of  h-,  arc  tabulated  in  column  7  of  Table  I.  In  the 
case  of  C0„,  for  T  <  500°R,  the  value  of  f^ (t)  used  in  computing  h1t  is 
differ*  .t  from  the  one  listed  in  colur.ui  8.  This  selection  wan  made  in 
order  to  reduce  the  error  in  the  approximation.  Values  of  the  enthalpy 
obtained  by  the  above  procedure  were  compared  with  values  in  the  Gas 
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Tables  at  various  temperatures .  The  error  vac  found  to  be  snail - less 

than  one  percent  in  most  of  the  eases,  (exceptions  were  IIg  at  lov  temper¬ 
atures,  vhere  the  error  is  about  2  percent,  and  CC2  in  the  neighborhood 
of  500°R,  v/here  the  error  is  found  to  be  4  percent).  These  error r,  be¬ 
come  negligible,  in  viev  of  the  large  values  of  the  heats  of  reaction. 


The  equilibrium  constant  K_  (T )  is  discussed  in  Fart  II,  vhere  the 
folloving  equation  is  given: 

AW1 


/  W'-  n, 

VT) '  r  7 


The  function  !•’  (t)  nay  be  obtained  from  tabulated  values  of  the 
equilibrium  constant,  K^,  given  in  terns  of  pi'cscures  (see  G.G.  Fenner, 
"Chemistry  Problems  in  Jet  Propulsion").  It  van  found  that  a  good  fit 
to  these  values  is  obtained  for  both  the  H„0-  and  the  C0o-  reactions 
vhen  F,  (t)  is  expressed  as  a  fifth  degree  polynomial.  In  the  range  of 

iw 

1200-3000°K,  the  error  in  the  values  of  1C  (t)  is  about  0 .1  percent 

7'k 

for  the  H20-rcaction,  and  about  3.0  percent  for  the  C02-rcaction .  'flic 
polynomials  ore  given  in  column  11  of  Table  I.  The  variable,  T  ,  shown 
there  is  a  function  of  the  dimensionless  temperature,  T,  and  is  given 
by  the  relation 

T'  -  uk  [im  - 120°] 


The  factor  1.30  is  used  to  convert  the  temperature  to  degrees  Kelvin, 
an  required.  It  vac  pointed  out  in  Part  II  that,  belov  approximately 
l800°R,  the  equilibrium  constant  becomes  large  enough  that  the  ratio 
k'/k  may  be  considered  zero.  It  is  unnecessary,  therefore,  to  provide 
n  curve -fit  for  temperatures  belov  this  limit . 

The  values  of  Ahf°,  the  heats  of  reaction  listed  in  column  6,  voro 
obtained  from  UBS  Circular  500  (Selected  Values  of  Chemical  Thermodynamic 
Properties,  February  195<- ) >  and  arc  given  at  absolute  zero  for  the 
component  in  the  gaseous  state.  The  values  of  B.  and  E . ,  given  in  column 

■X*  ^  * 

5  and  6,  arc  ponevhat  uncertain.  Because  of  the  assumption  that 


^c. 

These  values  were  suggested  by  F.  Boynton,  member  of  the  Applied 
Research  Group,  Convair  Astronautics 
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oxidation  of  H  ,  and  CO  proceed  according  to  a  single  r.tep  chemical 
reaction  (discussed  in  Section  II1-I32  of  Part  II )  it  bee  one  a  difficult 
to  deduce  values  of  and  from  the  meaner  experimental  data  avail¬ 
able  at  the  present  time.  For  this  reason  1?  has  born  grouped,  along 
with  other  reference  quantities,  into  a  single  parameter,  Kp  ,  an  shown 

lv 

in  the  above-mentioned  section.  The  solution  of  the  problem  ran  thus,  be 
carried  out  parametric  ally  for  various  values  of  K,  in  the  ranee  of 
variation  of  D  . 

B.  Reference  Values  and  Dimensionless  Parameters 


'flic  reference  values  used  in  non-di mensionali sin"  the  flow  parameters 
may  be  expressed  arbitrarily  in  terms  of  the  initial  values  in  the  jet 
and  the  air  stream.  One  nay  use  the  mean  values  for  the  velocities  and 
temperatures.  Thus, 


U 

o 


U  +  U 

e  a 


(27) 


For  a  reference  pressure, 
nozzle: 


T  +  T 


o 


(28) 


one  may  use  the  pressure  at  the  exit  of  the 


p  =  p  (Jbs/in2) 
o  e  ' 


(29) 


and,  for  a  reference  molecular  v/eight,  the  molecular  weight  of  the  jet 
mixture  at  the  nozzle  exit: 


(50) 


A  reference  density  is  obtained  from  the  equation  of  state, using  the 
reference  values  given  above.  Thus, 


P 


o 


T  R 

o 


(51) 


In  consistent  units,  this  equation  is  rewritten  as 

a  \i 

i|>4_  zn- 

1 

o 


po  13  hp  A 5  T  ^  ^ 


(31a) 
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The  reference  specific  heat  is  set  equal  to  the  specific  heat  of  the 
mixture  at  the  nozzle  exit.  The  specific  heat  of  each  species  at  the 
exit  is  given  by 


(32) 


and  the  reference  specific  heat  is  then 


c  =  c  =  /  (c. )  (c  ) 
po  pc  L  ie  pi  n 
i  L 

This  equation  can  be  modified  Doncvhat  by  using  the  relation 


(55) 


(Cj  - 


Vi 


Vi 


1  C  V  Wn 

4  Vi  ° 

Substituting  also  Eq.  (52),  Eq.  (55)  becomes 

%-vr  Lxi  fiQ 


(5»») 


o  i 


where  the  function  f ^ (Tc/Tq)  can  be  computed  from  column  B  of  T  bile  I. 

The  dimensionless  parameters  which  are  functions  of  tin  r>  fcrcncc 
values  given  above  will  be  rewritten  here  in  consistent  units.  Those  -’re, 


RT 


KP 


o  P  U  V,'  U 
00  rao  0 


or 


j. 

K.  =  52 .17  x  15li5  A5  — 


(55) 


4,  u- 


0  o 


and 


ir 

3  =  ° 


o  c  T 
no  o 


or,  in  consistent  units, 


U 


3  = 


o  32.17  x  778.26  c  Tq 


(56) 
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Tvo  other  dimensionless  parameters,  appearing  in  the  expression  for  the 
rate  of  production  of  species,  may  also  be  included  in  t'nir.  list.  The 
dimensionless  frequency  factor  is  giver.  by  [nee  Pcu't  II,  l’q .  (ph)] 

guyX>rK 
Uo  (xj 

x 

or,  in  consistent  unite,  expressing  B  in  (cn3)  /(inolc )2/(°R)  /sec,  and 
p  in  lbs/ ft3,  an  Given  by  Eq.  (31a), 

P2  B. 

-  5 .1520  x  10"1  — L  -*! 

K)  ° 

where  the  neccr.cary  convernion  factor  bar,  been  used  in  order  to  convert 

the  dennity  to  grm/cm3  an  required.  Finally,  the  dimensionless  activation 

*  /  \ 
energy  E^,  Given  by  Eq.  v25)  of  Port  II 


* 


K, 


in  written  in  the  following  form  when  E  in  Given  in  col/gr-molo, 

*  K 

T  in  °R: 

3  v 

1»0  * 

TTPol  tq 


and 


(3ft) 


C .  Input  Data 


In  addition  to  the  thermodynamic  data  Given  in  Table  I,  the  follow¬ 
ing  input  data  are  required  for  the  solution  of  the  problem.  The 
varioun  terms  arc  dieted  in  the  order  in  which  they  appear  in  the  computer 
program.  This  program  -trill  be  discussed  in  Section  III. 


Mc  Mach  number  of  the  Jet  at  x  =  0 

T\.  Ratio  of  specific  heats  for  the  jet  cide 

Ttpj  stagnation  temperature  of  the  jet  (°R) 

free  stream  or  flight  I  Inch  number 

M  Mach  number  on  the  air  side  at  ::  =  0 
a 

T  ratio  of  specific  heats  for  the  air  side 

e 

°  / 

T  free  stream  temperature  (°n) 

00 


lb 


pTJ 

dl 

k 

PTe 


0 

II 

(C1) 


X 


p 


ctognation  pressure  of  the  jet  (lbs/ in") 
lower  (air-side)  boundary  of  the  nixing  layer 
mesh  size  in  the  x-direction 
stagnation  pressure  of  the  air  stream  (lb/in2) 
reference  length,  nozzle  c::it  radius  (ft) 
turbulent  Prandtl  number 
turbulent  Levis  number 
Jet  spread  factor 

number  of  internal  grid  points  (olvnys  an  even  number) 

concentrations  (mass  fractions)  of  species  "i"  on  the 
r.ir  side,  x  =  0 

mole  fractions  of  species  "i"  on  the  Jet  side  at  =  0 

coefficients  appearing  in  the  polynomial  for 
li,  [Eq.  (9)1 

ing  in  the  polynomial  for 


a  ,a  ,a 

1'  2*  3*  4 


,P ,,P  ,(3  ,P  coefficients  appear! 
C0  2  3  4  K  [Eq.  (15)J 


The  above  coefficients,  as  veil  as  II  and  II  ,  are  obtained  from  the 
potential  flov  solution  of  the  Jet1  and  arc  discussed  in  Section  I.  The 
majority  of  the  data  arc  conditions  on  the  jet  or  the  air  stream  that 
muGt  be  specified  for  each  particular  ease .  L  and  k  may  be  selected 
arbitrarily .  A  discussion  on  the  value  of  DL  is  given  in  Section  III 
below. 


The  values  of  the  turbulent  Levis  and  Prandtl  numbers  arc  not  veil 
knovn.  Recent  c;:?erinental  data  on  recovery  factors  have  indicated  a 
value  of  Pr,  0.86,  and  this  value  nay  be  used  for  lack  of  other  data. 
In  any  ease,  valuations  in  the  Prandtl  number  are  not  ejected  to  be 
large.  Data  on  the  turbulent  Levis  number  is  also  lacking.  The  indica¬ 
tions  arc  that  itG  value  ranges  from  1.0  to  1.4.  It  vill  be  useful  to 
cany  out  the  solution  of  the  present  problem  for  various  values  of  Let 
in  thiG  range,  in  order  to  determine  its  effect  on  the  variation  of  the 
flov  parameters . 

The  Jet  spread  factor,  o,  required  in  the  t  ran  g  format  ion  of  the 
^-coordinate  is  not  veil  known.  Experimental  values  of  0,  and  its 
possible  variation  with  i'ach  number  arc  discussed  in  the  Appendix. 
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Finally,  it  should  be  pointed  out  that,  h,  the  ctep  Disc  in  the 
♦-direction,  does  not  appear  in  the  above  lict.  It  ic  computed  from 
1)^  and  II,  ao  discussed  in  Section  IV  of  Port  II  • 

III  .  COMPUTER  PROGRAM 

The  equation!)  derived  in  Part  II,  along  with  the  information  con- 

•x 

tained  in  Sections  I  and  II  of  this  report,  were  programmed,  using  the 
Fortran  II  system,  for  calculations  on  an  IBM-704  computer.  The  program 
consists  of  a  number  of  subroutines  which  are  employed  ns  required  by 
the  procedure.  These  subroutines  contain  numerous  descriptive  comments 
and  arc,  in  general,  self-explanatory.  They  vill  be  listed  here,  and 
discussed  briefly,  for  case  of  reference. 


Subroutine 


INPUT:  provides  for  storage  of  input  data,  listed  in  Section  II-C 

CONST:  used  to  calculate  initial  conditions,  given  in  Section  I-A, 

and  the  reference  values  of  Section  II-B 

PREKE:  executes  preliminary  operations,  such  as  assigning  the 

initial  and  boundary  conditions 


GEJJV: 


DVAL: 


FSUBI,FSIP, 
GIT,  and  IIBI: 

KOMKN,SPSCI, 

EIIEEG: 


contains  the  calculations  for  functions  of  x,  such  as 
Mach  numbers,  pressures,  pressure  gradients,  and 
the  variable  F,,(v.J  [which  contains  the  reaction  term] 

K.  Iv 

used  for  calculating  boundaty  values 

used  to  calculate  thermodynamic  functions  f. (t),  f*(T), 
hn  ,  and  G±  (t  )  listed  in  Table  I.  1  1 

used  to  calculate  all  coefficients  of  momentum,  species, 
and  energy  equations . 


SOLVE:  used  to  solve  the  above  equations  simultaneously 

IIICRE:  incrementing  subroutine  used  to  obtain  the  solution  of 

the  equations  at  each  successive  step  in  the 
x-direction 


TRAIiSF:  program  used  for  carrying  out  the  transformation  of  the 

♦-coordinate  to  the  physical  coordinate  y 


The  program  vac  prepared  by  D.  Witte  and  P.  Doherty,  of  the  Digital 
Computer  laboratory,  Convair,  Astronautics,  and  was  assigned  the 
number  595-42042 . 
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An  additional  subroutine,  denoted  as  subroutine  "MEA1J5"  on  the  first 
pace  of  the  program,  is  planned  for  the  future .  This  subroutine  \;ill 
make  it  possible  to  conputc  the  coefficients  of  the  equations  nt  a  raid- 
point  of  a  mesh,  rather  than  at  the  preceding  point,  thus  imp  roving  on 
the  accuracy  of  the  numerical  solution. 

The  number  of  internal  grid  points  is  limited  to  a  maximum  of  lO, 
because  of  the  limited  memory  of  the  704  computer.  The  value  of  II  =  l3 
vould  require  the  solution  of  l44  equations,  and  the  coefficient  matrix 
of  these  equations  would  occupy  (f)  percent  of  the  memory  cells,  leaving 
approximately  11000  cells  for  storage  of  other  quantities,  and  for  the 
planned  subroutine  discussed  in  the  preceding  paragraph .  It  is  possible 
by  means  rf  on  external  tape,  to  increase  the  capacity  of  the  memory  so 
that  values  of  II  larger  than  l8  may  be  accomodated.  However,  the  com¬ 
putation  tine  will  then  be  increased  from  4  minutes,  to  approximately 
20  minutes,  for  each  step.  It  nay  become  necessary,  at  some  future  date, 
to  use  the  external  tape  for  at  least  a  fev  runs  in  order  to  be  able  to 
evaluate  mesh  sice  errors. 

The  program  in  its  present,  form  allows  for  only  one  value  of  the 
lower  boundary,  I3L,  to  be  specified  at  the  start  of  the  solution.  The 
larger  the  value  of  D_,  the  further  doimstrcam  can  the  solution  be  car- 
ried.  However,  for  large  values  of  13^,  the  step  sice  in  the  V  direction 
?s  large,  and  it  becomes  difficult  to  obtain  accurate  profiles  at  the 
initial  stages  of  the  mixing  layer.  It  is  ncccssaiy,  then,  to  start 
with  a  snail  value  of  BT,  obtain  the  solution  for  the  first  few  steps, 

i-l 

and,  then,  restart  the  program  vith  a  larger  value,  in  order  to  carry 
the  solution  further  downstream. 

To  remove  this  difficulty,  it  is  planned  to  incorporate  into  the 
program  a  nev  subroutine  which  will  cause  the  value  of  BT  to  increase 
automatically,  as  the  solution  progresses.  The  increase  vill  take  place 
whenever  the  value  of  the  temperature  at  the  boundary  of  any  given  step 
reaches  the  usual  limit  of  99  percent  of  its  asymptotic  value . 
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IV. 


TYPICAL  CASE  AlfD  PREUMUIARY  RESULTS 


Ac  a  typical  problem,  vre  will  concidcr  the  cace  of  a  two- dimen Dionol 
Jet  of  Mach  number  Mjet®  3.185  exhausting  into  an  air-stream  of  Mach 
number  1^,  *  5*70  at  an  altitude  of  150, 000  feet.  Table  II,  belov,  gives 
the  boundary  of  the  Jet  and  the  Mach  numbers  and  total  pressure  at  each 
Dido  of  the  boundary,  ao  shown  in  Figure  la.  The  coordinates  X/re  and 
r/rc  are  ao  shown  in  Figure  lb.  The  angle  is  the  angle  that  the  tangent 
to  the  boundary  makes  wi-n  the  axiG  of  symmetry  of  the  nozzle,  and  is  not 
used  here.  The  dimcnBionlcso  coordinate,  x  lor  x  ),  used  in  this  work 
was  related  to  the  coordinates  of  Table  II  by  Eq.  (8)  aB  discussed  in 
Section  I-B .  The  following  simple  relation  wbb  found  to  hold  to  within 
*1.5  percent: 

—  =  0 .345  x 
r 


The  Mach  numb  ere  M.  and  M  were  exprccscd  qb  fourth  degree  polynomial  0 

J  ® 

in  x#  Thus 

H.  -  a  +  a  x  +  a  xs  +  a  x3  +  a  x4 

J  0  1  2  3  4 


where  the  values  of  the  coefficients  ore 


ao  -  4.86946  ' 
ax  *  .28303 

«2  *»  -.64624  x  10’ 1 
a3  =  .85049  x  lo"2 
«  - .42125  x  10"3 


The  otrean  side  Mach  number  is  given  by 

M  aP  +  0,x  +  0  :;2  +  0  X3  +  0  ::4 

S  CO  1  2  3  4 

with  the  following  values  for  the  coefficients 

\o  =  1 .36321 

0X  «  .29241 

»  -  .68381  x  10" 1 
03  *  .91005  x  10"2 

P4  *  -  .43438  X  10’3 
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TABLE  II* 


DATA  OBTAINED  FROM  THE  INVISCID  JET  FLOW  SOLUTION1  FOR  JET 
MACH  NUMBER  M  =  3.185  AND  ATTITUDE  150,000  ft 


CHAR 

x/r 

e 

r/r 
'  e 

Degrees 

J1 

O.OCO 

1.000 

33-53 

J2 

0.090 

1.059 

33*35 

SO? 

0.134 

1.088 

33.26 

J3 

0.179 

1.118 

33.32 

SU3 

0.191 

1.126 

33.51 

SD3 

0.205 

1.135 

33.30 

S03 

0.270 

1.178 

33.23 

j4 

0.358 

1.235 

33.17 

J5 

0.537 

1.352 

33.00 

S07 

O.817 

1.533 

32.807 

S09 

1.096 

1.712 

32.51 

SD11 

1.357 

1.877 

32.21 

SD16 

2.024 

2.29? 

31.67 

j8 

2.651 

2.674 

30.95 

J9 

3.815 

3.186 

30.17 

J10 

4.373 

3.680 

29.65 

Jll 

5-223 

4.160 

29.15 

J012 

6.081 

4.633 

28.66 

SD4b 

6.529 

4.876 

28.36 

SU48 

6.820 

5.032 

28.17 

PoJ 

lb/ in2 , ABS 

M 

s 

^os 

lb/ in2, ABS 

4.872 

492.0 

1.965 

3.1 

4.896 

492.0 

1.091 

3.1 

4.907 

492.0 

2.002 

3.1 

4.928 

492.0 

2.024 

3.1 

4.946 

492.0 

2.042 

3.1 

4.936 

492.0 

2.03? 

3.1 

4.955 

492.0 

2.051 

3.1 

4.981 

492.0 

2.079 

3.1 

5.028 

4o2.0 

2.126 

3.1 

5.096 

492.0 

2.196 

3.1 

5.147 

492.0 

2.248 

3-1 

5.187 

492.0 

2.288 

3.1 

5.285 

492.0 

2.387 

3.1 

5-340 

492.0 

2.441 

3.1 

5.412 

492.0 

2.513 

5.1 

N“\ 

CO 

J- 

lA 

492.0 

2.584 

3-1 

5.538 

492.0 

2.638 

3.1 

5.591 

492.0 

2.691 

3.1 

5.614 

492.0 

2.713 

3.1 

5.628 

492  -0 

2.726 

3.1 

This  Table  was  furnished  by  Dr.  J.M.  Rowyer,  member  of 
the  Applied  Research  Group,  Convair,  Astronautics.  It 
was  computed  by  the  method  of  characteristics. 
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The  error  in  the  above  approximations  is  Iccc  than  0.5?>  of  the  tabulated 
value o. 


In  addition  to  the  above  coefficients,  and  the  thermodynamic  data 
riven  in  Table  I,  the  following  input  data  were  used:  (the  cyribols  ore 
defined  in  Section  II-c) .  Mc  ■  4.872,  TJf  *  1.225,  *  6300% 

t'v,  -  5.700,  Mn  «  1.965,  rE  -  1.40,  Tw  -  505%  pTJ  -  492.0  lbo/in2. 


PTo  a  5.075  lbs/ in2,  L  -  1.91  ft,  Prt  -  O.QO,  %  =  1 
W  -  Ik,  Dl  -  -  1.50  and  -3.0,  k  -  0.125. 


.20, 


56.0, 


The  concentrations  of  epeeies  at  x  =  0,  and  at  the  upper  and  lower 
boundaries  arc 


Species 


(Ci> 


h2o  0  0 .566 

H2  0  0 .1420 

o2  0 .253  0 .0072 

co2  0  0 .1372 

CO  0  0 .3476 

He  0 .767  0 


The  value  of  BT  =  -1.50  was  used  to  compute  the  first  five  steps  in 

the  development  of  the  mixing  layer  (x  =  0.625).  A  run  was  then  made 

vith  Bl=  >3.0,  and  the  analysis  was  carried  to  x*  1.25,  i.e.,  1-1/4 

nozzle  exit  radii.  The  results  are  plotted  in  Figures  2  through  6,  in 

terms  of  both  the  ♦  and  the  y  coordinate.  The  computed  reference  values 

are  U  *=  7436.6  ft/sec  and  T  ■  1926. 5  *R.  The  velocity  and  concentration 
o  o 

profiles,  Figures  2,  3,  and  6  have  the  expected  trend,  and  are  not  af¬ 
fected  by  the  value  of  B^.  It  is  noted  in  Figures  4  and  5  that  the  tem¬ 
perature  profiles  exhibit  a  "bulge"  near  the  x-axis.  The  peak  thus  ob¬ 
tained  is  lower  by  about  %  for  the  larger  value  of  BL  (and  therefore  h). 
It  is  believed  that  this  error  is  the  result  of  the  larger  n,  sh  sizo. 

* 

The  values  of  X,  listed  above  arc  based  on  the  results  obtained  by 
Hcu  and  Boynton.2 
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FIGURE  4b  DIMENSIONLESS  TEMPERATURE  PROFILES 

C 
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FIGURE  60  CONCENTRATION  PROFILES  DURING  THE  INITIAL 
STAGES  OF  DEVELOPMENT  OF  THE  MIXING  LAYER,  x*  =  0  25 


•0.07  -0.05  -003  -001  0  001  003  005 


FIGURE  6b  CONCENTRATION  PROFILES  DURING  THE  INITIAL 
STAGES  OF  DEVELOPMENT  OF  THE  MIXING  LAYER,  x*  =  0.50 


2k 


1 


y* 

FIGURE  6c  CONCENTRATION  PROFILES  DURING  THE  INITIAL 
STAGES  OF  DEVELOPMENT  OF  THE  MIXING  LAYER,  x*  *  075 


y* 

FIGURE  6d  CONCENTRATION  PROFILES  DURING  THE  INIT'AL 
STAGES  OF  DEVELOPMENT  OF  THE  MIXING  LAYER,  x*  =  1.125 
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In  order  to  determine  whether  the  bulge  In  the  temperature  profiles 
develops  into  a  sharp  peak,  and  also  to  obtain  the  thickness  of  the  mix¬ 
ing  layer  further  downstream,  the  solution  was  carried  to  x  =  5«'>  (5  nozzle 
radii)  using  B  = -12.0.  A  plot  of  the  computed  thickness  is  given  in 
Figure  7*  To  reduce  the  computation  time,  the  value  of  k,  the  step  size 
in  the  x-direction  was  doubled.  This  increase  introduced  some  error  in 
the  peak  values  of  the  temperature  profiles  (not  reported  here).  The  net 
increase  between  the  peak  value  and  the  asymptotic  value  at  the  lower 
boundary  was  found  to  be  l4$  at  x  =  5*0  as  compared  to  12.5$  at  x  =0.5. 

It  appears  that  the  cooling  of  the  gases  at  the  boundary,  as  a  result  of 
the  expansion,  and  the  low  value  of  the  pressure  in  the  mixing  layer 
(p  =  0.17  lb/in^  at  x=5*0)  prevent  the  bulge  in  the  temperature  profile 
from  developing  into  a  sharp  peak. 
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FIGURE  7  COMPUTED  THICKNESS  OF  THE  MIXING  LAYER 


NOMENCLATURE 


Kjj  dimensionless  parameter  defined  by  Eq.  (37) 

K 

L  reference  length 

M  Mach  nunber 

II  nunber  of  internal  grid  points 

p  pressure  ( (linen Gionle sc  when  not  subscripted,  otherwise  has 
the  units  of  lbs/ in2) 

p  diracnGionloss  pressure  gradient 

r  y  coordinate  of  the  jet  boundary  as  sho'.m  in  Figure  lb. 

R  universal  gas  constant 

R  gas  constant  for  the  mixture 

T  temperature  (dimensionless  or  in  °R  aG  noted) 

U  velocity  at  the  boundaries  (ft/sec) 

W  molecular  v, 'eight  of  the  mixture 
x  dimenGionless  coordinate 

X  x  coordinate  of  the  jet  boundary  aG  shown  in  Figure  lb. 

X.^  nolo  fraction  of  specieG  i 

GubGcriuts 

. .  .  <r,  ,,  , 


a 

e 

j 

o 


T, 


T 


denotes  initial  conditions  (::  =  0)  for  the  air-stream  Gide 
•initial  conditions  (x  =  0)  for  the  jet  Gide 
jet-side  conditions  for  x  >  0 
reference  values 

conditions  on  the  air- stream  side  for  x  >  0 
total  or  stagnation  conditions  on  the  jet  side 
total  conditions  on  the  air- stream  side 
free  stream  conditions 
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C*4  coefficients  appearing  in  the  polynomial  for  Mj 

P  .  coefficients  appearing  in  the  polynomial  for  M 

CO  '  G 

T..  l'atio  of  Gpccific  heats  for  the  Jet  Gide 

X  ratio  of  specific  heats  for  the  air  side 

o  Jet  spread  fact'' - 
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APPENDIX 


The  value  of  o  is  to  be  obtained  from  experiments  on  the  rate  of 
spread  of  supersonic  Jets.  At  the  present  time,  the  available  data  are 
limited  to  values  of  M  <  1.60,  and  the  results  of  various  investigators 
vary  widely,  as  seen  in  Figure  8.  In  this  figure,  the  quantity  a refers 
to  the  value  of  o  for  an  incompressible  Jet  (m  ■*  0).  The  value  of 

=  12,  has  been  established  by  numerous  experiments.  For  M  >  2  values 
of  the  ratio  a/aQ  obtained  from  the  experimental  data  of  Anderson  and 
Johns,  (jet  Propulsion,  Vol.  25,  No.  1,  1955)  on  Jets  from  solid  propel¬ 
lant  rocket  motors  have  been  used,  for  lack  of  other  data.  For  M  =  1.40 
and  1.84  the  data  on  hot  air  Jets  in  this  reference  is  in  fair  agreement 
with  the  results  of  other  investigators.  A  smooth  curve  has  been  drawn 
through  the  various  points  to  describe  the  possible  variation  of  a/aQ 
with  M.  Obviously  this  curve  is  only  tentative  and  should  be  revised 
when  more  information  becomes  available 
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FIGURE  8  EXPERIMENTAL  DATA  ON  THE  VALUE  OF  <r 


